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Entanglement in Spin Waves with Binary Dipolar Atoms in a 2D Optical Lattice 
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We study the area-dependent entropy and two-site entanglement in the spin waves with binary 
interacting dipolar Bose-Einstein condensates in a 2D optical lattice. The two species condensates 
at each site behave like an ensemble of spin-half particles with the interaction to its nearest neigh- 
bors and next nearest neighbors. We show that the ground state of a Bose-Einstein condensate 
lattice with nearest-neighbor and next-nearest-neighbor interactions exhibits spin wave and can be 
mapped into a harmonic lattice. The critical behavior of the area-dependent entropy and two-site 
entanglement is examined in the vicinity of a quantum phase transition. 

PACS numbers: 03.75.Gg, 03.75.Lm, 03.75.Mn 



Since Bose-Einstein condensates (BEC's) of alkali 
gases were first observed in a magnetic trap the ex- 
perimental study of BEC's has grown rapidly. In particu- 
lar, multi-component BEC's have been realized using the 
different hyperfine states of 87 Rb [2] and spinor conden- 
sates using different Zeeman states of sodium F = 1 0- 
In addition, interspecies interactions give rise to interest- 
ing phenomena such as phase-separation £| . Besides this, 
the quantum phase transition (QPT) 5] from a superfluid 
phase to a Mott insulator phase has been observed us- 
ing atoms in optical lattices @- BEC's in optical lattices 
are also being used to study aspects of quantum entan- 
glement [jj in many-body systems because of their con- 
trollability 8J. Indeed, evidence of multi-particle entan- 
glement has been seen in atoms with two internal states 
prepared in Mott insulator in optical lattices 

In this Letter, we study the entanglement in the spin 
waves of a two-component BEC .2\ in a 2D optical lat- 
tice |9| in which the condensates at each site behave like 
a "spin magnet" and interact with each other to its near- 
est neighbors and next nearest neighbors. Spin waves 
can be obtained in this spin lattice in the low excitation 
regime. Moreover, this "long-ranged" interaction can be 
produced by dipole-dipole interactions, and with this in 
mind we should note that Bose-Einstein condensation of 
chromium, with a high magnetic dipole-dipole interac- 
tion, has recently been realized [ljj. This dipole-dipole 
interaction strength can be tuned by ma gne tic fields, or 
by engineering the geometry of the trap \\\\ . Such con- 
trollable interspecies, short- and long-ranged interactions 
are important in producing the multi-particle entangled 
states. 

The study of the entanglement measures leads to a 
novel perspective on the structure of ground state and 
their critical behavior. We shall show that this "BEC 
lattice" is equivalent to a set of harmonic oscillators in 
the low-excitation regime. This means that the well- 
known analysis of a quantum harmonic lattice can be 
used to investigate entanglement in this many-body sys- 
tem. We shall further see that the entropic measure of 
entanglement is useful to understand our system's prop- 
erties of area dependence which used in some parts of 



quantum field theory [12]. In the case of the Klein- 
Gordon (KG) field, the entropy of the field is obtained by 
tracing over the variables outside a region under consid- 
eration and the entropy found to be directly proportional 
to the boundary area of inside region. This was shown 
by Bombelli et al. and Srednicki by considering a free, 
massless and scalar KG field, and which is also equivalent 
to the ground state of a coupled of harmonic oscillators 

m 

An analysis of the entropy and entanglement in a ID 
harmonic lattice system was from the viewpoint of quan- 
tum information theory was given in the reference |14| . 
along with the generalization to the 2D and 3D cases. 
The area-dependence of the entropy fl5| and the quan- 
tum correlations [16j in a harmonic lattice including some 
aspects of the critical behavior were also studied. Never- 
theless, the general relationship between area-dependent 
entropy and the nature of QPT is not yet fully under- 
stood [15| . We will examine an aspect of this relationship 
and provide some new insights. 

Two-site entanglement is also a most useful quantity 
with which to examine the nature of quantum correla- 
tions in a lattice system. For the case of an infinite spin- 
chain, with nearest neighbor interaction, the critical be- 
havior of two-site entanglement was given in references 
fl?| . It is clear from these studies that two-site entangle- 
ment shows the non-local nature of a many-body system 
in a QPT. 

In this Letter, we consider a two-component BEC 
trapped in a square optical lattice with MxM sites as 
shown in Fig. The atoms involved have two different 
magnetic states with tunable magnetic dipolar interac- 
tion strength 0HI3. We take this dipole-dipole 
interaction strength between atoms to be proportional 

to 1/r 3 Enm, 

where r is the atoms between two sites. 
Moreover, we assume that these two internal magnetic 
states of atoms can be coupled by applying an external 
magnetic field. We adopt the single-mode approxima- 
tion ^3 to the condensates trapped deeply in each site 
in which the condensates can be described by the local- 
ized mode functions associated with the potential wells. 
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FIG. 1: The configuration of the 2D lattice with the nearest- 
neighbor interactions gi and Q2 in the horizontal and vertical 
directions respectively, and with the next-nearest-neighbor in- 
teraction 2~ 3 ^ 2 g2- 



The Hamiltonian of the system is then given by 



H 



Ha- 



(1) 
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<i,J> 



+^(4a 3 + a\ ai ) + <>;;';/,>, + 6t ftj)j (2 ) 
Y^uKAh + bU)^, (3) 



where and hi are the annihilation operators of compo- 
nents A and B in the i-th site. The interaction parame- 
ters k^(k^), n l ab and u)i are respectively the intraspecies 
interaction, the interspecies interaction and the coupling 
rate between the two magnetic states in the i-th site. 
The hopping and the dipole-dipole interaction strengths 
between the i-th site and the j-th site arc denoted by 
rt l J(ClV) and g l i — g\ or g 2 for the component A(B) re- 
spectively. We consider the atoms in each site to interact 
with the nearest neighbors and the next-nearest neigh- 
bors as shown in Fig. ^ 

For the purposes of discussion, we consider the interac- 
tion parameters and the number of atoms to be the same 
for each site, i.e., lo % — ui, k — (n l a + n\ — «^ & )/4 and 
N l = N. In the Mott-insulator limit, £1<Ck, the tunnel- 
ing between the sites is negligible and hence the trapped 
two-component condensates in each site can be regarded 
as an ensemble of pseudo spin-half particles. For this 
case, we can write the Hamiltonian in terms of angular 
momentum operators in the following form: 



H x = 



4kJ! 



<ij> 



1 ji jj 



(4) 



Here J l x = (a\a z - b\bi)/2, J l y = (ojftj - 0*6^/2*, J\ = 

{a\bi + a,ib\)/2 and iVj = aja; + b\bi is the total number 
of atoms in the i-th site. 

We assume a sufficiently large external magnetic field 
is applied to the system such that u> <C ft, g and in the 
order of kN. We can then apply the Holstein-Primakoff 
transformation (HPT) |2(j to map the spin operators into 

the harmonic oscillators, where J]_ = c\\l N — c\a 



\J N — c\aci and J\ — (c\ci — N/2). In the low degree of 

excitation regime, (cjci)/A<Cl, the Hamiltonian can be 
approximated thus: 



(5) 



<i,3> 



This effective Hamiltonian H 2 is the zero order approxi- 
mation of this exact Hamiltonian Hi and gives a descrip- 
tion of spin waves Q. In the thermodynamic limit, this 
approximation becomes exactly equivalent to the system 
H\. The exact numerical solution of this multi-spin sys- 
tem even for just a few sites is extremely difficult. This 
analytical although approximate solution is therefore ex- 
tremely valuable and provides important insight into the 
physics of this multi-spin system. Indeed, the approxi- 
mation turns out to be valid in the two-site case even if 
the number of atoms N is in the order of several hundred 
|2l| . We immediately see from this that the ground state 
and low-lying eigenstates of our system will behave like 
a set of coupled harmonic oscillators. We can see this 
by representing the operators in the position-momentum 
space: qi = (cj + c 4 )/\/2w and pi = iy/oj(c\ - Cj)/\/2. 
Then, the Hamiltonian can be written explicitly in terms 
of a set of coupled harmonic oscillators: 



H' 2 



(6) 



<hJ> 



J'_ 



Here, is the potential matrix that expresses the inter- 
action strength between the oscillators qi and qj . In the 
infinite lattice limit, the spectrum of the Hamiltonian can 
now be found by using a two-dimensional Fourier trans- 
form. 

The onset of a second-order QPT can now be deter- 
mined by simply finding when the energy gap A, i.e., the 
energy difference between the first excited and ground 
states of the system, vanishes [f|. This model is now 
directly solvable, and so the spin wave modes and the 
energy gap A can be evaluated explicitly. This means 
that we can also determine the critical coupling g 2 as a 
function of g\. We find the following results 

[ (4k - 2 91 + u/N)/(2 - V2), g 2 < V2gi 
92 = \ (4k + u/N)/2, _ gi = V2gi 

[ (4k + 2 5 i + uj/N)/(2 + V2), g 2 > V2gi 

(7) 

Thus, the phase diagram within this Gaussian approxi- 
mation can be determined directly from Eq. (7f . Here, 
we consider the case when interaction strengths g\ and 
g 2 are equal, i.e., g\ = g 2 = g. The critical cou- 
pling g c is (lu + AkN)/N(A — \ 2). As g approaches to 
g c , the lowest excitation energy gap is given by A = 
[ivN(4:-^2)] 1 ^ 2 \g-g c \ 1 ^ 2 . If g is greater than the critical 
value g c , then the excitation gap A becomes complex, the 
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low-lying excitations approximation breaks down, and 
the Gaussian approximation thus fails. This implies the 
pseudo spin-half particles in each site will be "polarized" 
to some extent which can be viewed as "quantum magne- 
tization" of the spin-half atoms |9( . We intend to examine 
this issue in our future work. 
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FIG. 2: The entropy E L i as a function of L is shown, for 
w = 500k;, AT = 1000 and the lattice size MxM are 6400. 
For the infinite (finite) case, the different coupling strengths 
g is approximate to g c — 1.74028 ... up to 11 decimal place in 
unit of k, 1.5 k, 1.25 n are shown, where are denoted by solid 
line (empty circle), dashed line (empty square) and dotted 
line (empty up-triangle) respectively. 

We can now examine the entanglement properties of 
the ground state of our system. As the ground state is 
a Gaussian state, its 2M x 2M density matrix can be 
completely determined from the second order moments 
(XiXj + XjXi) — 2(Xi)(Xj), where Xj are the quadra- 
ture variables qi{pi) and % = 1, . . . , M. The density ma- 
trix of the system p can be in terms of the covariances 
(XiXj) only since (Xj) = in the ground state. 

We investigate the area-dependent properties of the 
entropy of the bipartite entanglement of two sets of har- 
monic oscillators. We consider the MxM square har- 
monic lattice is bisected into a L x L square lattice and 
the remaining M 2 — L 2 lattice as "inside" and "outside" 
parts respectively, say 1 and 2. Moreover, we consider a 
LxL "inside" square lattice which is located at the cen- 
ter of the whole lattice whereas the rest of oscillators are 
called the "outside" part. The bipartite entanglement 
can be determined by the von-Neumann entropy E L i (p) 
of the reduced density matrix which is given by 

E L2 (p) = -Tr( Pl log 2(0l ), (8) 

where pi = Tr 2 (p) is the reduced density matrix of p by 
tracing out the "outside" subsystem 2. The reduced den- 
sity matrix p\ can be obtained by including the position 
and momentum correlations in the set of modes 1 only. 
It is given by 



where Q and P are the position and momentum covari- 
ance matrices with the matrix elements Qij = (qiqj) and 
Pij = (piPj) of the subsystem 1 respectively. The entropy 
measure of this bipartite entanglement between two re- 
gions E L 2 is then found to be [l!| 




(10) 

where the symplectic eigenvalues are the square root of 
the eigenvalues of QP or PQ. It is important to note that 
all symplectic eigenvalues must be greater than or equal 
to 1. Also, we follow the definition of the bipartite entan- 
glement in the reference 15] in which they count the dou- 
bly degenerate symplectic eigenvalues only once to show 
the entanglement area law being valid in a finite-ranged 
interaction harmonic system with arbitrary dimension. 

We now proceed to examine the relationship of the 
entropy Ei? and the boundary length L. In fact, by using 
Fannes' inequality [22^. we can show that this bipartite 
entanglement entropy of the many-body Hamiltonian Hi 
and that of harmonic system Hi are very close for a finite 
block L [2j|. Hence, the harmonic system is a good model 
to study the area law of this exact many-body system. In 
Fig. |21 the entropy E L i against the boundary length L 
with different couplings g are shown. The entropy Ejj. is 
directly proportional to the length L in the non-critical 
regime. Moreover, in the case of the coupling g closing to 
the critical coupling g c is also shown to satisfy the area- 
law. Therefore, the entropy of the bipartite entanglement 
is still area-bounded in the vicinity of the QPT. For the 
finite-size case, we consider the total number of sites are 
MxM = 6400 in the system as shown in Fig. [21 From 
our numerical results, they are clearly shown to be a good 
agreement with the infinite case and also satisfy the area- 
law. 




s 



FIG. 3: The derivative dC,i/dg as a function of g (in units of 
At) is shown with the same parameters to the previous figure. 
The infinite case (solid line), M = 41 (dashed line), 31 (dash- 
dot line) and 21 (dotted line) are shown respectively. 
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Next, we study the entanglement of atoms trapped in 
two different sites. The problem of two-site entanglement 
is to determine the entanglement between any two har- 
monic oscillators in an ensemble of harmonic oscillators. 
Indeed, it is analogous to determining two-mode entan- 
glement in quantum optics The criterion of the in- 
separability of two-mode system [24| can be applied to 
the pure and the mixed states and hence it can be used 
to evaluate two-site entanglement. The two-site reduced 
density matrix p(*> 5 'J can be constructed from the correla- 
tion functions of these two sites i and j only. By applying 
several local transformations, the reduced density matrix 
of a symmetric Gaussian state can be written as |2 ll 

/ n c \ 

»(y) - [ n o -c j , m 

p c n I ' [ ' 

\ -c n ) 

where n = 2{(qf) (p?)) 1 /* = 2«< ? f >(pf ))V2 and c = 

^(^(lil^iPiPj)) 1 ^ 2 ■ This two-site entanglement param- 
eter can be defined as C\i-j\ = n — c |2l|. If Cli-jl i s 
below one, then it is said to be entangled. Moreover, this 
parameter can be used to evaluate the amount of entan- 
glement of formation (EOF) j2f|. This amount of EOF 
is a function of the parameter C\i-j\ f° r < C|<— j'l < 1 

mm. 

In the infinite limit, we can treat our system as a sym- 
metric Gaussian state due to the preservation of transla- 
tional symmetry. Therefore, we are able to evaluate the 
EOF through the two-mode entanglement parameter in 
this case. In the finite-size case, we investigate the oscil- 
lators locating at the center in the square lattice and the 
total number of oscillators are odd. This Gaussian state 
is nearly symmetric if the system size is large enough. 
We study the two-site entanglement parameter £i of two 



adjacent sites. It decreases as the strength g increases 
as shown in the inset in Fig. [3] but the minimum is not 
formed at the critical point. This means that the two- 
site entanglement can be controlled by the strength g. In 
addition, we found that it is only the nearest neighbors 
that are entangled. The quantum entanglement is thus 
very short-ranged indeed. In fact, this feature coincides 
with the intuition of the area-dependent entropy. It is be- 
cause the oscillators in the "inside" region entangle with 
the oscillators in the boundary. The entropy is thus at 
most proportional to the boundary L. 

In order to examine how the ground state changes close 
to a QPT [17| . we numerically examine the first deriva- 
tive of the two-site entanglement d(i jdg of the nearest 
neighbor. In Fig. [31 d(i/dg diverges as closing to the 
QPT in the infinite case whereas it also increases greatly 
to the critical point g c in the finite-size cases. Moreover, 
the larger size cases increase with a larger rate with re- 
spect to g. This result obviously shows the finite-size 
scaling in the critical regime. 

In conclusion, we have investigated the entanglement 
properties in the spin waves of the two-component con- 
densates in a 2D square lattice with the finite-ranged 
interaction closing to a QPT. We show that the entropy 
still satisfy the area-law and the scaling of two-site en- 
tanglement is examined in this finite-ranged interaction 
system. The measurement of this Gaussian state here 
is envisaged by an experimental technique that has been 
devised recently to measure the correlation functions be- 
tween atoms |2g. Moreover, the physical realization of 
harmonic chain may lead to the applications of quantum 
information science p7|. 
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